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Inelastic neutron scattering studies have been carried out on selected phonons and magnetic ex¬ 
citations of a crystal of superconducting (SC) CaioPt 4 As 8 (Fei_xPtccAs)io with the onset transition 
temperature ~ 33 K to investigate the role that orbital fluctuations play in the Cooper pair¬ 
ing. The spectral weight of the magnetic excitations, at Q = Qm (magnetic P points) is 

suppressed (enhanced) in the relatively low (high) uj region. The maximum of the enhancement is 
located at cj = cdp ~ 18 meV at temperature T = 3 K corresponding to the ~ 6.3. 

This large value is rather favorable to the orbital-fluctuation mechanism of the superconductivity 
with the so-called 5'++ symmetry. In the phonon measurements, we observed slight softening of the 
in-plane transverse acoustic mode corresponding to the elastic constant Cqq . This softening starts at 
T well above the superconducting Tc, as T decreases. An anomalously large increase in the phonon 
spectral weight of in-plane optical modes was observed in the range of 35 < cj < 40 meV with de¬ 
creasing T from far above Tc. Because this cj region mainly corresponds to the in-plane vibrations 
of Fe atoms, the finding presents information on the coupling between the orbital fluctuation of the 
Fe 3d electrons and lattice system, useful for studying possible roles of the orbital fluctuation in the 
pairing mechanism and appearance of the so-called nematic phase. 


I. INTRODUCTION 

In the study of physics of Fe-based superconductors^ 
it is important to understand what roles the orbital fluc¬ 
tuations plays in realizing the high superconducting (SC) 
transition temperature Tc. Of course, it is widely known 
that various characteristics of the systems, two dimen¬ 
sional layers of strongly correlated electrons and super¬ 
conductivity in the neighbourhood of the magnetically 
ordered phase, suggest that the spin-fluctuation exchange 
is relevant to the SC pairing^i^ and that the so-called S± 
symmetry of the order parameter A is realized with its 
sign reversal between the Fermi surfaces around F and 
M points in the reciprocal space {Q space). However, 
the situation is not simple, because the possibility of the 
5'++ symmetry without the sign reversal was pointed 
out based on the observed small rates of the Tc sup¬ 
pression by non-magnetic impurities much smaller than 
those expected for the S± symmetry^-— It has been also 
shown theoretically that the orbital fluctuation mecha¬ 
nism can explain the high-Tc superconductivity with the 
5'++ symmetry as well as various other experimental re¬ 
sults which were simply regarded as evidences for the 
S± symmetry A”— Among these works, the primary roles 
of the orbital fluctuation enhanced by the Aslamazov- 
Larkin type vertex correction due to the spin fluctuation 
is emphasized. 

Experimentally, the phase diagram of 
Ba(Fei_a;Coa;) 2 As 2 has attracted much attention^iii 
because we can find many physical characteristics which 


are common to many other Fe-based systems. One 
example is the existence of the tetragonal-orthorhombic 
structural transition at Ts higher than the antiferro¬ 
magnetic transition temperature Tn, implying that the 
orbital degrees of freedom should be taken into consid¬ 
eration besides the spin degrees of freedom. Another 
example is that even in the macroscopically tetragonal 
phase, the breakdown of the four fold symmetry be¬ 
comes appreciable with decreasing T, in static quantities 
such as the electrical resistivityi^ and energy splitting 
of the electronic band formed by the 3dyz and 3dzx 
orbitalsi^i2£ at around the characteristic temperature 
(T*) significantly higher than Ts (we call T* nematic 
temperature here.). In arguing these findings, it seems 
to be inevitable to have roles of orbital fluctuation in 
mind, as stressed by Zhang et al^ for example. Now, 
to answer the question “ 5 pm or orbital? ” has become 
primarily important for the understanding of the physics 
of strongly correlated electrons entangled with spins, 
orbitals and probably lattice systems. We have carried 
out neutron scattering studies on magnetic excitations 
and phonons of a CaioPt 4 As 8 (Fei_a^Pta^As)io crystal 
with the onset SC transition temperature ^ 33 

K to find effects of the orbital fluctuation on these 
dynamical behaviors of Fe-based systems. In the phonon 
measurements, we focused on particular phonons near 
the F and M points in its Q space of the pseudo tetrag¬ 
onal cell. These reciprocal points were chosen, because 
the orbital fluctuation is expected to be strongJ^d^ 
Parts of the results of studies have been published in 
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FIG. 1. (a) Unit cell of CaioPt 4 As 8 (Fei-xPtccAs)io is shown 
schematically, (b) Diamagnetic signals of a part of the sam¬ 
ple. The closed and open circles show the result in field and 
zero-field cooling conditions, respectively, where the field was 
applied parallel to the c direction. 


Refs, iiii 


II. EXPERIMENTS AND CALCULATIONS 


CaioPt 4 As 8 (Fei_a;Pta^As)io was discovered by Kakiya 
et alM It has the stacking of-(Fei_a;Pta;As)-Ca-PtAs 2 - 
Ca- layer units as shown schematically in Fig. [IJa). 
Within the a-b planes, the so-called ^/bao x ^/bao unit 
cell is formed, uq (^3.903 A) being the lattice parameter 
of pseudo tetragonal cell of the plane s^^i^^ having two Fe 
atoms in an FeAs plane. The PtAs 2 planes were found, 
by NMR 1 /Ti ^ to be essentially insulating and nonmag¬ 
netic, and although the crystal structure is rather com¬ 
plicated, we can study the physical properties of the FeAs 
conducting planes of CaioPt 4 As 8 (Fei_a^Pta^As)io, mainly 
considering the electronic structure similar to other Fe- 
based superconductors ^^^22 

A single crystal, which weighs ^ 10 g, was prepared 
from the starting mixture of CaAs : FeAs : Pt : Fe with 
the molar ratios 1.00 : 1.00 : 0.51 : 0.09 (Ca : Fe : 
Pt : As = 1.84 : 2.00 : 0.94 : 3.67). The susceptibility 
measurement on a part of the crystal with a SQUID mag¬ 
netometer exhibits a significant shielding diamagnetism 
below ^ 33 K as shown in Fig. [D^b). The dis¬ 

tance between the FeAs layers (or lattice parameter c) of 
the pseudo tetragonal unit cell estimated from the OOL 
reflections (Fig. [2]) is ^ 10.51 A. In the same figure, we 
can clearly distinguish the epitaxially grown second phase 
CaioPt 3 As 8 (Fei_a;Pta^As)io with the FeAs-layer distance 
c of ~ 10.33 A. This difference between the c values of 
these two kinds of phases is commonly found in many pa¬ 
pers published previously. Judging from the very small 
fraction of the Bragg reflection intensities of the second 
phase, we can find that the second phase does not affect 
results of the present studies on their dynamical proper¬ 
ties. As shown in the Appendix, the x value is estimated 
to be 0.05-0.06, which is smaller than the value reported 
previouslyThe rapid growth of the shielding dia¬ 
magnetism with decreasing T through the SC 

transition occurs in the major part of the sample just 
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FIG. 2. X-ray diffraction intensity of OOL reflections of a part 
of the crystal used for the neutron scattering measurements. 
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below or at least, within the T range of a few 

degree below We will return this point later in 

relation to the magnetic excitation data. In the studies 
of the NMR Knight shiftand 1 /Ti, X-ray structural anal¬ 
yses, magnetic susceptibility and electrical resistivity, no 
evidence for structural and antiferromagnetic transitions 
has been observedwhich is contrasted to the case 
of CaioPt3As8(Fei_3;Pt3;As)io4ii^i2^ 

Data of the spin excitations and low-energy phonons 
were collected with both the thermal (2T1) and cold 
triple-axis-spectrometers (4F2) at the neutron reactor 
ORPHEE of Laboratoire Leon Brillouin (LLB), Erance, 
respectively, where incident and scattered beams were 
focused by pyrolytic graphite 002 monochromater and 
analyzer, respectively. The collimation conditions were 
fully open. Neutron data of optical phonons were col¬ 
lected with the thermal triple-axis spectrometer (HB-3) 
at HEIR of Oak Ridge National Laboratory, USA. The 
collimation condition was 48'-40'-40'-70'. In all cases, py¬ 
rolytic graphite filters were placed before the analyzer to 
eliminate the higher order reflections. 

The phonon energies have been calculated for 
CaioPt 4 As 8 (EeAs)io by the first-principle calculation 
package, Vienna Ab initio Simulation Package (VASP)^ 
where the structure of CaioPt 4 As 8 (EeAs)io optimized 
before the phonon calculation was used (The crystal pa¬ 
rameters thus obtained are in the previous report 4^; 
space group P4/n). The calculated phonon results are 
compared with the experimentally observed ones. 


III. RESULTS AND DISCUSSION 
A. Magnetic scattering 

Eirst, we briefly show results of magnetic 
scattering obtained for the crystal sample of 
CaioPt4As8(Eei_a:Pta;As)io. The typical Q- 

scan profiles of the magnetic scattering intensity 
^mag(Q,^) = (n + 1) Xmag(Q 7 ^) obtained at constant-cj 
with the scattered neutron energy Ef = 14.6 meV are 
shown in the left panel of Eig. [3l n being the Bose factor, 
where the directions of the longitudinal and transverse 
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FIG. 3. (Left) Q scan profiles of the magnetic excitation spectra measured along L- and T-directions (inset), respectively, 
through the M point Q = (0.5, 0.5, 0) at T = 38 K. (Right) The difference of the intensities of the magnetic excitations at the 
M point Q = (0.5, 0.5, 2) between T = 3 K (SC) and 38 K (normal) (solid circles), and between 24 K (SC) and 38 K (normal) 
(open circles). The arrow indicates the energy of the peak position of the difference. 


scans through the M point (1/2, 1/2, 0) (L- and T-scans, 
respectively) are also shown in the inset. The anisotropy 
of the profile width is clearly observed as reported 
commonly for various other Fe-based systemsThe 
difference of the magnetic scattering intensities between 
38 K (> Tc) and two temperatures, 24 K and 3 K (< Tc) 
observed at another M point Q = (1/2, 1/2, 2) are in 
the right panel of Fig. [3l From the figure, we can see the 
following. (1) With decreasing T through Tc, or with in¬ 
creasing SC order parameter A, the magnetic scattering 
intensity in the low-energy region is suppressed and it is 
enhanced in the high-energy region, and the maximum 
of the enhancement is found at cj = cjp ^ 18 meV, which 
corresponds to c^p/ZcbT®^®®^ ^ 6.3. (2) The spectral 

weight was found to be zero in the low cj-region (at 
least cj < 7 meV at 3 K), indicating that the SC state 
is established above 24 K all over the whole volume, 
as stated in section IT The peak of the enhancement, 
is usually expected for SC order parameters A with 
the sign reversal, and called “the resonance peak”. For 
Fe-based systems, as Maier and Scalapino^ pointed 
out, it appears, if the spin-fluctuation mechanism is 
relevant to the superconductivity, at cjp smaller than 
the value of |Ar| -b |Am|, Ap and Am being the order 
parameters at the Fermi surfaces around the F and M 
points, respectively. However, it was also pointed out 
in Refs. [I and 0 that even for the *5++ symmetry of 
A, such a peak can be expected at cjp > |Ar| + |Am|, 
if we consider the electron dissipation effect. On this 
point, we cite Ref. liO, which shows that the observed 
relation uj^/k^Tc ^ 4.3 holds for various systems, and 
that larger values of uj^/k^Tc have been found for 
CaioPtnAs8(Fei_a;Pta^As)io (n = 3 and 4) systems, 
including the present data c^p/ZcbT®^®®^ ^6.3. Because 
it is not easy to experimentally determine the A values 
accurately, we cannot distinguish which one of the S± 


and 5++, the relation co’p/ZcbT®^®®^ 4.3 itself supports. 
Then, the relation c^p/ZcbT®^®®^ ^ 6.3 observed here, 
which is favorable to the 5'++ symmetry, is rather 
encouraging. (Note that cJp/ZcbT®^®®^ observed here is 
regarded as its minimum estimation in the sense that 
we used the onset value instead of Tc.) Then, 

we proceeded the phonon measurements, the results of 
which are described in the next section. 


B. Phonons 

Phonon measurements have been carried out for two 
kinds of modes with main characters of the in-plane 
transverse acoustic (TA) phonon near the F point and 
in-plane optical modes at an M point mainly polar¬ 
ized within the FeAs plane, because the orbital fluc¬ 
tuation is theoretically expected to be strong at these 
pointsThe former phonon corresponds to the elas¬ 
tic constant Cqq , which is well-known to exhibit the soft¬ 
ening in Ba(Fei_a^Coa^) 2 As 2 as T decreases toward the 
tetragonal-orthorhombic transition point Ts^^i^ In the 
present case, the corresponding TA phonon was mea¬ 
sured at several Q = (2,0,0) + q and T points with 
q = (0,g',0). The observed cj-scan profiles are shown 
in Fig. lUat two fixed g^-values of 0.05 and 0.25, for ex¬ 
ample. They were obtained with the neutron energy- 
loss condition and with the final energy being fixed at 
Ef = 5.57 meV and Ef = 8.00 meV, respectively. Be¬ 
fore the measurement at g = 0.05, the Bragg point (2, 
0, 0) was adjusted at each T point to ensure an accurate 
value of q. The phonon energies were determined by fit¬ 
ting the Gaussian line shape as shown in Fig. IH where 
the intensity is divided by the Bose population factor 
n + 1. The results are shown in Fig. [5] in the form of 
ujq{T)/ujq{T = 250 K) for q = 0.05, 0.15 and 0.25. (Cal- 
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4F2/LLB, CalO-4-8, Q = (2, q, 0), TA phonon 
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FIG. 4. The TA phonon profiles obtained by scanning the 
transfer energy uj at Q = (2,0,0)+ with two q vectors of 
(0, 0.05, 0) and (0, 0.25, 0) are shown. The solid curves are 
fitted lines. The horizontal lines are their baseline, shifted 
arbitrarily for each temperature from 5 K to 260 K. 


culations using the convolution of the resolution function 
gave essentially the same results.) In the figure, we can 
see the softening tendency of the phonon energy oX q = 
0.05 as compared to those at the higher q points with 
decreasing T, that is, the energy at q = 0.05 shifts to 
the lower-energy side, while those at q = 0.25 have the 
opposite tendency. The relative magnitude of the soften¬ 
ing fraction extrapolated to (g, T)^0 is consistent with 
the data of the elastic constant Cm reported for the opti¬ 
mal or slightly overdoped Ba(Fei_a,Coa,) 2 As 2 ^ which is 
also SC and exhibit none of the magnetic and structural 
transitions. (Here, note that uj oc \fCm-) 

The dispersion curves of this in-plane TA branch was 
compared with those obtained by the first principle cal¬ 
culation for the pseudo tetragonal cell, where the calcu¬ 
lated and experimental results were found to agree rea¬ 
sonably well , even though we used the crystal structure 
of CaioPt 4 As 8 (FeAs)io optimized by the calculation. 

We next move to the study of the in-plane optical 
phonons. The phonon intensity was measured at the 
fixed scattering vector Q = (2.5, 0.5, 0) by scanning 
the energy transfer cj, where the final neutron energy 
was fixed Ef = 14.7 meV. Note that this spectrome¬ 
ter setting is to see the phonon modes mainly having 


4F2/LLB, Ej^ = 5.57 meV, Q = (2, g, 0), TA phonon 
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FIG. 5. Temperature dependences of the phonon frequencies 
of the TA branch measured at Q = (2,0,0)+ g are shown 
for q = (0,g, 0),(g = 0.05,0.15, and 0.25) in the form of 
ujq{T)/ujqiT — 250 K)—T plot. The curves in the figure are 
guides for the eye. 


in-plane polarizations at the M point. Then, after sub¬ 
tracting the cj-independent background intensity (B.G.) 
determined in the high energy region from the raw data 
we have obtained the phonon spectral weight 
Xph(QM,^) = [/ph(QM,^) - B.G.]/(n -h 1), and the re¬ 
sults are shown in Fig. [6] against uj for various T values. 
We can find significant T dependence of Xpj^((3M,^) in 
the region of uj around 37 meV. Note that no spectral 
weight exists above 40 meV, as shown below by the cal¬ 
culation. The data for cj > 45 meV are just due to the 
excess B.G. at Q = (2.5, 0.5, 0). 

To study the origin of this strong T dependence, 
we first calculated, by the first principle calculation^ 
the partial phonon density of states g^'^\Q^uj) = 

Y^- e^j\Q) S{uj - uJjiQ)) with e^'^(Q) being the i (x, 

or z) component of the polarization vector of the j-th 
mode (or we often specify the atomic element of the j-th 
mode, too, by only this j). The results are shown for the 
i-component of Fe and As motions in Figs. [7|^a) and[7|^b), 
respectively, where we find that the peak positions of the 
total density of states seem to be slightly shifted to the 
lower side from the observed peak positions at 300 K. 
We can also see that all modes at the M point are in the 
region less than 40 meV and the peak at the largest en¬ 
ergy originates mainly from the in-plane vibrations of Fe 
and As atoms. Although the observed phonon intensity 
(Fig. [6]) cannot be compared directly with the calculated 
data shown in Fig. [71 partly because the former has the 
factor I (Q) • Q P instead of e^j\Q) in the latter, we 
can, at least, say that the observed highest-energy peak 
of g^'^\Q,uj) corresponds the in-plane motions of FeAs 
layers. 

From this view point, we plot in Fig. [8] the observed 
phonon intensity integrated over the uj region of 35-40 
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FIG. 6. Phonon spectral weight at the M point Q = (2.5, 
0.5, 0) is shown for several T points. These spectra are cor¬ 
rected by Bose factor after subtracting the oj-independent 
background determined in the high energy region. 
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FIG. 7. Galculated phonon density of state defined as de¬ 
scribed in the text is shown at M point for (top) Fe atoms 
and (bottom) in-plane As atoms. Each component of the vi¬ 
brational direction can be distinguished by the colors of the 
line. The total density of states represents the sum of the 
ones with i — x, y, and z. Note that the red and green lines 
almost overlap with each other. 


meV against T, together with those in the regions of 10- 
16 and 19-23 meV, where the difference among the T 
dependences of XphiQ'^^^) explicitly seen. The data 
at Q = (2.5, 0.25, 0) are also shown by crosses taken at 
5 K and room temperature. 

Here, we add the following point. First, the anoma¬ 
lous T dependence of the phonon intensity cannot be ex¬ 


plained by considering the so-called “false peaks”, which 
may appear accidentally due to the spectrometer condi¬ 
tion (see Ref. for details). Second, effects of the Debye- 
Waller factor on the T dependence of phonon intensities 
is almost independent of the phonon energy uj^ at a fixed 
Q point, which is also confirmed by the first principle cal¬ 
culation in the present case at Qm = (2.5, 0.5,0). The 
calculation has also confirmed that the Debye-Waller ef¬ 
fect is much smaller than the observation at around oo ^ 
37 meV, but seems to be consistent with those in the 
other uj regions. 

Now, we have found the gradual but rather significant 
increase of the phonon spectral weight of the in-plane op¬ 
tical modes. It starts at ^250 K with decreasing T, and is 
reminiscent of the similar T dependence of various static 
quantities observed in Fe-based systems in the region be¬ 
tween the structural transition (Ts) and nematic (T*) 
temperatures (Note that Ts = 0 in the present 

sample.) In this T region, the breakdown of the four-fold 
symmetry, in particular, the splitting of the ^dyz and 
^dzx orbitals have been foundi^i^S suggesting the contri¬ 
bution of the strong orbital fluctuation to stabilize the 
local lattice distortion. Then, the observed softening of 
the in-plane TA phonon and the unusual T dependence of 
the in-plane optical phonons should contain information 
how the lattice system couples to the orbitals (and/or 
spins). In turn, it can be used to distinguish which one 
of “spins or orbitals” is primarily important for the real¬ 
ization of the nematic phase and superconductivity. 

It is worth noting that if phonons primarily contribute 
to the pairing mechanism, unusual T dependence of 
Xph(Q’^) expected as the superconductivity occurs at 
the phonon energy cj ^ 2A^ which is much smaller than 
the energy region of the observed unusual T dependence 
(around ^ 37 meV). This excludes the possibility of the 
primary relevance of the phonons to the SC mechanism. 
(Note that the orbital fluctuation mechanism itself does 
not require the primary role of phonons.) Finally, the 
strong T dependence of XphiQ^^) seem to be 

confined within the Q region very close to the M point, 
as shown by the data at (5 = (2.5, 0.25, 0), which remains 
as a future issue to be clarified by detailed studies. 


IV. SUMMARY 

We have carried out inelastic neutron scattering mea¬ 
surements of phonons as well as magnetic excitations of a 
single crystal of CaioPt 4 As 8 (Fei_a:Pta;As)io. The maxi¬ 
mum of the x((^g^g(QM, -enhancement by the supercon¬ 
ductivity is observed ai uj = ujp ^ 18 meV, which corre¬ 
sponds to co’p/ZcbT®^®®^ ^ 6.3, rather favorable to the 5'++ 
symmetry of the order parameter expected by the orbital 
fluctuation mechanism. In the measurements of the in¬ 
plane TA phonons, the softening at the F point has been 
observed, suggesting that it is a common phenomenon of 
Fe-based superconductors. In the measurement of the in¬ 
plane optical phonons at the M-point Q = (2.5, 0.5, 0), 
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FIG. 8. T dependences of spectral weights integrated over 
the regions 10-16, 19-23 and 35-40 meV are shown by open 
circles, open squares, and closed circles, respectively. Cross 
symbols show results of 35-40 meV region taken at Q = (2.50, 
0.25, 0). Data are scaled by the values at T = 300 K for each 
Q point. The curves are the guides for the eye. 

the strong T dependence of the phonon spectral weight 
Xph(QM,^) was found. These phonon behaviors become 
appreciable with decreasing T far above Tc. The results 
obtained here present information to understand what 
roles the orbital fluctuation plays in realizing the super¬ 
conductivity and also nematic behavior found in various 
Fe-based systems. 
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Appendix A: Pt content in the FeAs planes 


In our previous reports, the x value of the present sam¬ 
ple was estimated to be ^ 0.2 by analyzing the line shapes 
of ^^As NQR and NMR-centerline spectra^^ 

However, considering the data reported in Refs. ii,® 
and in, we have checked the x value again, and found 
that, by using the value x = 0.06, the NQR and NMR 
data can be explained better than the previous case, as 
shown in Fig. A[T1 Because it is consistent not only with 
the existing reportsbut also with the used start¬ 
ing composition, we revise the value. By this revision, 
no essential qualitative changes of the results in the pre¬ 
vious reports emerge except the ma gnit udes of the NMR 
Knight shift shown in Fig. [3] of Ref, b?!. 
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FIG. A 1 . The NQR profiles at 40 K (top) and angular (0) 
dependence of the NMR profiles (bottom) taken at 100 K 
with the applied magnetic field H = 9.4013 T are shown with 
the fitted curves (solid lines) obtained for x = 0.06, 0 being 
the angle between H and the c axis of the pseudo tetragonal 
cell. The N value in AsN indicates the Pt number doped 
to the four nearest neighbor Fe sites, and 2c and 8g in the 
parentheses indicate the two crystallographically distinct AsO 
sites (space group: P4ln). In the calculation, the ratios of 
the signal intensities of As 0 ( 2 c), As 0 ( 8 ^) and Asl sites were 
approximated as 0.16 : 0.62 : 0.22 by considering the binom¬ 
inal distribution of Pt impurities doped to Fe sites of FeAs 
planes.— (Note that we did not consider the difference be¬ 
tween the two distinct sites of Asl( 2 c) and Asl( 85 ') , because 
it does not result in meaningful changes.) 
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